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SIMPLE MODELS FOR THE SHUTTLE REMOTE MANIPULATOR SYSTEM 


Introduction 

One of the fundamental objectives of the NASA Space Shuttle Mission 
is to develop the capability to launch and retrieve payloads from near-earth 
orbit using the reusable shuttle orbiter vehicle. One of the more demanding 
of the early space shuttle missions will be the deployment of the 
Interium Upper Stage (IUS) which is currently under development for the 
Department of Defense. The primary function of the IUS will be to place 
a single large spacecraft or multiple smaller spacecraft in earth orbits 
which differ from the shuttle orbit or to place planetary spacecraft on 
earth- escape trajectories. The IUS-vehicle will be approximately sixty feet 
long and cylindrical in shape with a diameter of about fourteen feet. 

The overall weight of the IUS can be as much as 60,000 lb, which will be 
about one third of the overall shuttle-orbiter weight. The IUS will rely 
on a three-axis stabilized propulsion and avionic system for trajectory and 
attitude control. As a consequence, it is critical that the IUS spacecraft 
module be deployed in a stable attitude for check-out and launch. In the 
event of malfunction during the deployment stage, it will be necessary 
that the IUS module be recovered and returned to the shuttle payload bay. 
Hence, the RMS must be able to accomodate each of these operations. 

The requirement that these manuevers be conducted while the shuttle 
and IUS are in close proximity without disturbing the attitude of the IUS 
vehicle, lead to several critical operational requirements. One of the 
primary problems associated with this close proximity operation arises 
from the IUS interactions with the shuttle RGS thruster plume. Any 
translational motion required by the shuttle for station-keeping and/or 





2 


rendezvous purposes must be performed using the 900 lb RCS thrusters. 

Since the projected area of the 60 feet long by 14 foot diameter IUS 
module will be on the order of 8A0 square feet, a substantial force c.an 
be generated by the pressure of the RCS jet plume. Since, in general, 
the center of pressure will not coincide with the center of mass, this 
interaction will lead to significant moments on the IUS. Tumbling and/or 
severe disturbances in the IUS attitude may occur. If these moments 
fire Introduced during the retrieve stage, severe forces and moments may be 
introduced in the remote manipulation system RMS/arm. Such moments may 
lead to structural failure of the RMS arm. 

In addition, a careful study of the effects of the shuttle RCS thruster 
actuation during any of the close proximity maneuvers must be made. The 

primary questions of concern are those related to the effects of the moments 

■» 

on the IUS module. the forces generated on the RMS arm, control strategies 
and limitations on the RMS arm movement to prevent damage and, finally, 
thruster firing limitations to prevent unacceptable moments on the IUS. 

The basic RMS design calls for two manipulator arms and the supporting 
equipment for the operation of the arms. One manipulator arm will be mounted 
on the port side of the payload bay on all missions. The second arm will 
be installed on the starboard side of the payload bay and will be removable, 
if it is not required during a given mission. The RMS will be composed of 
the following major components: (.a) two manipulator arms, (b) the RMS dis- 
plays and control, and (c) the manipulator control interface unit. These 
components will be supported by the associated RMS software as well as the 
RMS ground control equipment. As with all space rated systems, weight and 
power requirement must be kept to a minimum. This leads to severe design 
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requirements which can be met if all pertinent factors are included in the 
design process. 

While the state-of-art in the design and utilization of fixed-based 
remote manipulator systems has achieved a comparatively high degree of 
sophistication through medical and industrial applications, the specific 
applications of this technology to the shuttle remote manipulator system is 
not a direct extension of the current practice. The fundamental difference 
between the shuttle problem and those previously treated lies in the fact 
that the shuttle will be operating in an orbiting environment. As a 
consequence, the shuttle and the associated payload, which is either to be 
deployed or retrieved, will be interacting in an essentially zero-gravity 
environment. In this situation, unbalanced forces developed due to the mutual 
shuttle-payload interactions will cause rotational motions about the center 
of mass of both the shuttle and the payload. Since the RMS arm will be required 
either to move to accomodate these relative motions or to counter the motion 
by elastic deformation, severe design restraints are placed on the operation 
and control of both the shuttle and the RMS in the payload deployment/retrieval 
mode. Pre-launch study of the IUS operation can be accomplished only through 
numerically simulated studies. The simulation models will be complex since 
a number of factors will contribute to the phenomena under study. However, 
the models and the associated computer simulation program will serve as basic 
mission analysis tools. Consequently, it is mandatory that NASA develop 
a comprehensive understanding of the characteristics, accuracy, and limitations 
of the computer program simulations. 

Computer programs which give reasonable complete descriptions of the 
motion are described in Refs. 2 and 3. The computing effort required to 
simulate motions of the arm with these models is quite severe and hence 
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case studies cannot be made in real time. In addition, there is significant 
computation costs associated with exercising computer programs based on complete 
mathematical models. As a consequence, simple models which contain the 
essential characteristics of the dynamic motion are of interest for pre-mission 
design studies. The models would be useful for performing a number of case 
studies rapidly with minimal computer costs. The conclusions reached after 
study with such models would then lead to a limited number of special cases 
which would be studied in detail by the more precise mathematical models. 

With this objective in mind, this investigation is aimed at establishing a 
series of simple models which can be used to study the forces and moments 
which occur due to the RCS jet plume firings during a deployment or retrieval 
of an IUS type payload. The models considered in this investigation are 
primarily planar in nature. The extension to three-dimensional motion is 
straightforward and will be addressed in subsequent studies. In this study 
primary attention is given to the roles the payload play in determining the 
overall moments on the arm. 



Approach 


In order to gain insight into the overall dynamics, the planar model 
of the space shuttle/RMS/payload system shown in Fig-. 1.1 is considered. 

As stated previously, the primary objective of the analysis will be to 
characterize the forces and moments at the shoulder of the RMS (point a) as 
a function of the thruster force, T. The analysis will be developed in 
several stages. In the first stage, the shuttle/RMS/payload combination 
will be assumed to be totally rigid. After considering the equations of 
motion for this model, a model with a linear spring at the shoulder, a, 
will be developed. The model will then be extended to include linear springs 
at both a and b. In each analysis, force and loading moments at a will 
be identified. Following a complete study of the single joint model, similar 
studies will be made of a two-joint model in which the arm is separated into 
two segments connected by linear spring at the midpoint C. 



Figure I. 1: Shuttle/RMS/Payload Model 

The analysis will be developed in several steps. In the first step, the 
Shuttle/RMS/Payload combination will be assumed to be totally rigid. Next, 
a model with a linear spring at A will be developed. Then a model with linear 
springs at A and B will be considered. In each analysis, forces and bending 
moments at A will be sought. 
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I. RIGID BODY MODEL 

In this model, the shuttle, arm, and payload are assumed to be a 
single rigid body. The geometry is shown below in Figure 1.2. 



Figure 1.2: Rigid Body Model 



or, in terms of the motion of the center of mass, C, of the shuttle. Analyses 
based on G and C will be presented for comparison and to lay the groundwork 
for the more complex models to follow. 
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Motion Referenced to System c.m, ,G : 

The free-body diagram for the system is shown in Figure 1.3. 



Figure 1.3: Rigid Body Model Coordinates 

The overall goal of the analysis is, as stated earlier, to derive 
forces and bending moments at the shoulder, point A. The pertinent equations 
are: 



(m + m ) 
s p 


(m + m) 
s p 



^G ~ *GQ X ^ I G 01 


( 1 . 2 ) 


Since the selection of a thruster gives the components of T, and the specified 
shuttle-payload geometry gives r , equations (1.2) form three equations in 
the three unknowns a , a , and a. 

Cjr b 

X Z 

The geometry needed to determine r^ is simple since r G> r c , and r c are 


known 
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\ 


r Q 


r c + r cq 



(1.3) 


To find the forces and bending moments at A induced by the thruster firing, 
a free body diagram of the shuttle alone is necessary (see Figure 1.4). 



Figure 1.4: Shuttle Free Body Diagram 


For the shuttle. 


ZF = A 

+ T 

= m a 

X X 

X 

s c 

ZF = A 

+ T 

= m a 

z z 

z 

s c 

< 

II 

+ x 

CA 

X A + 


x 


(1.4) 


CQ 


yy„ 


a 


The unknowns are A , A , M , a. n , and a . From kinematics, the acceleration of 

X Z n L L 

X z 

C is related to the acceleration of G by 


S c = 5 G + “ x ( “ x *00* + 5 x r GC 


(1.5) 
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For the two-dimensional model, this vector equation results in two 
scalar equations. One new unknown is introduced, u). However, 0 ) is related 
to a by a = dai/dt. Thus, a value is needed only for U) q . 

In scalar form, the equations necessary for this analysis are as follows. 
Let the origin of the chosen coordinate system be located at the initial 
location of the system center of mass and be translating along with it. When 
a thruster is fired, of course, the center of mass, G, will move away from 
the origin of the coordinates. Thus, 


r_ = Oi + Ok 

° (I.! 

V„ = Oi + Ok 

Vjr 

O 

The analysis is much simplified if the thrust firing is represented by an 
impulsive change. Equations (1.2) become, considering equations (1.6) 


V G 

1 

(m + m ) 

</Tdt) 

X 

s p 


V G 

1 

(m + m ) 

</T dt) 
z 

z 

s p 


to = 

O' 

u,° 

r-||w C 

+ 

o 

3 

X /Tdt) 


In the third of equations (1.7), 



but initially r_ = 0. Also 

VJ 


(1.7) 
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Also 


and 


r q ’ r c + r CQ 

? q = (x c + x cq )I + (z c + z cq )E 
" V + Z Q 5 


T « T i + T k 
x z 


X G Q * f “ < 2 q T x ' x q T z* 


CX.8) 


Thus, if a constant thrust over a time of 0.1 is assumed, equations (1.7) become 


0. IT 


G (m + m ) 
x s p 

0.1T 


G (m + m ) 
z s p 

. Vx - Vz 
“ = “o + io t: ■" 

Lt 


(1.9) 


Equations (1.9) determine the velocity of the system center of mass 
and the angular velocity at the end of any given impulsive firing of the 
thrusters. This information will be used to determine steady-state loads. For 
dynamic loads, as the thruster is fired it is necessary to return to Equations 


( 1 . 2 ). 

Equations (1.2) can be rewritten as 

T 

x 


G (m + m ) 
x p s 

T 

z 

a G “ (m + m ) 
z p s 

a = z q T x ~ Vz 

*0 


(1.10) 
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These three equations give a„ , a_ and a explicitly. Considering oo as a 

b b 
X Z 

known parameter for the analysis which can be varied to determine its effects, 
the acceleration of the shuttle center of mass, C, can be formed using 
equation (1.5). For this purpose, let 
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Equations (I. 10), (I. 11), (1.12), and (1.13) form the computation sequence 
for A , A and M in the rigid body case for any thruster. The parameters 

X Z A 

which must be specified are: 


Choice of Thruster: 

(T , T , 
x’ z 

X CQ* Z CQ 

Payload Location and Mass: 

(m , x » 
P c 

Z C’ I yyp 

Shuttle Parameters: 

(I , m 

) 


yys 

s 


I G is given by geometry, I , and I yyp 


Overall System Parameters: 
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II. SINGLE JOINT MODEL 

The forces and torques acting on the shuttle are as shown below in 
Figure II. 1. 



where the quantities marked through are zero because of the assumption of planar 
motion. The xyz coordinate system is centered at the shuttle center of mass, C, 



14 


and rotates with the shuttle. The force A is the force exerted by the RMS 
on the shuttle at the shoulder and the force T is the force exerted on the 
shuttle by the firing of a thruster. In the current model, any thruster can 
be selected to exert the force T. The torque about C produced by firing the 
thruster is M . In this model, the torque exerted by the RMS on the shuttle 
at A will have only one component (in the j direction) . 

The angular velocity of the shuttle at any time, u), can be written in 
component form as 


0) 

s 


i + q j + r k 
s s 


For the two-dimensional model, this reduces to 


“ ■ v 


The force and moment equations for the shuttle in body-fixed coordinates 
are (for the two-dimensional model) 

(11. 1) 

(11. 2) 


ZF = T 

+ A = m (u + q 

X X 

X X n : 

ZF = T 

• 

+ A = m (w - q 

z z 

z x n 

ZM„ « M. 

+ r x a + M_ 

G A 

CA x T 

M. 

+ z„ . A - x„,A 

A 

CA x CA z 

I 

• 

a 

yys ‘s 


(II. 3) 


where m is the mass of the shuttle and I is the mass moment of inertia 
s yys 


of the shuttle about the y axis (through C) . 


From equations (II. 1) , (II. 2) , and (II. 3), we get 
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T + A 

* , x x , . 

U '-V +— m (I1 - 4) 

s 

T + A 

* - v + — ~ <«•« 

s 

% = I {M A + Z CA A X “ X CA A z + ^ (II. 6) 

yys 

Equations (II. 4), (II. 5), and (II. 6) would give the shuttle c.g. center 

of mass velocity and pitch rate versus time if they could be intetrated. However, 

A , A , and M. are unknown, 
x z A 

In order to determine A , A , and M. , it is necessary to consider the 

X Z il 

payload/RMS as a second rigid body and look at its motion. Let D be the center 
of mass of the payload/RMS combination. The shut tie /payload /RMS is shown in Figure 
II. 1 below, while the payload/RMS free-body diagram is shown in Figure II. 3. 



z 


Figure II. 2: Definition of <J) 
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-A 


Figure II. 3: Arm-Payload Geometry 


The equations of motion for the payload/RMS are 


EF - -A 

= 151 a TX 

(II. 7) 

X X 

P D 

r X 

EF = -A 

= m a^ 

(II. 8) 

y z 

P D z 

“O'-*! 

+ J DA x (“A) - I q 

P 

(II. 9) 


where m is the mass of the payload/RMS combination and I is its mass 

p yy p 

moment of inertia about an axis parallel to the y-axis through point D. 

Note that in equations (II. 7) and (II. 8) * a^ and a p are expressed in 

x z 

terms of a rotating coordinate system which is rotating at the angular rate 
of the shuttle (i.e., ^ coor( ^ s = = q s< j). 

From Figure (II. 2), note that 

r^ = d cos (j> i - d sin <j> k (II. 10) 

where <p is the angle between the positive x-axis and the line AD, and d is 

the distance from A to D. Thus, r CD is given by 

r = r + r 
CD r CA AD 

= ( X CA + d cos cf>)I + (z CA - d sin 4>)k 


(II. 11) 
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Ss 

The quantity d can be allowed to vary during relative motion of the shuttle 
and payload. The present analysis will allow this and it will be assumed 
that this motion is known over the period of interest and will be specified 
over the period of interest by three parameters, d o> d Q , and d Q . The motion 
will be modeled as 

d(t) = d +dt+d t 2 /2 (11.12) 

o o o 

The quantities d Q , d Q , and d Q will be input parameters in the computer program. 

It is now necessary to develop expressions for a^ and a^ . In general, 

x y 

i D = 5 c + i D +fix(nx r CD ) + 2Q x V D + H x ? CD (11.13) 
rel rel 

It is most convenient to develop each term of equation (11.13) separately and 
then combine them to form a^. In this manner, errors can be more easily 
avoided . 

First, a is known from equations (II. 1) and (II. 2), i.e., 

V# 

a r = (u + q w) i + (w - q u) k (11.14) 

L s s 

The quantity a^ is the acceleration of point D relative to point C disregarding 
rel 

coordinate rotation. Thus, a^ is the second derivative of r CD taken with 

rel 

unit vectors fixed. (It is convenient to do this in two steps as the first 

derivative will be identified as and will be used in a later term) . 

rel 
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Equation (II. 11) gives r^. Thus, 


\L = (+ d cost}) -d<£ sintj))i + (-d sint}) -d<|) cos<j>)k (11.15) 


since x_. and z„. are constant. Next, 
CA CA 


a * (+ d cos<}) -2dp sim}> — d<f> sin<}> -dc}> 2 cost}>)i 

rel 


• • • « — 

+ (-d sin<|> — 2d^> cos(}> -dcj) cost}) + d4> sin<}))k 


The quantity ft x (ft x r^) is given by 


Thus, 


q g j x [q g j x {(* CA + d cost}))! + (z CA - d sintfi)k) } ] 


S x (fi x ^ cd ) = (- q 2 - q 2 d cos<}>)i 


+ ( “ q s Z CA + q s d Sin * )E 


The term 2 ft x V_ is given by 
rel 


2ft x V = 2q j x { (+ d cost}) ~dp> sint())i 

*J m S 

rel 


+ (-d sin<{> — d<f> cost}))k} 


Thus, 


2ft x V n = (-2q d sint}> -2q dip cost}))i 
U ~ s s 

rel 

+ (-2q d cost}) + 2q dt}) sin<}>)k 
s s 


Finally, the term ft x r is given by 

Q x r CD = q g J X ^( X CA + d cos<}>)i + (z CA - d sint}>)k} 
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Thus, 


R x J CD * ( VCA ' q s d Sin ' t ’ )I + ( ' q s X CA ' q s d cos, W E 


(11.19) 


Combination of equations (11.14), (11.15), (11.18), and (11.19) results in 


1 

a^ = u + q g w + d cos<$ ~ 2dcf> sin (j> 
x 

- dtp sincj) - dip 2 cos <p - q 2 x 

S uA 

2 # 

- q d coscp - 2q d sin<j> 

s s 

~ 2q g d| cos (p + q g Z CA 

- q g d sincj) 


( 11 . 20 ) 


• • • 

a^ = w - q g u - d sin cp - 2d(p cos<j) 

z 

- dcj) coscj) + dcj) sin <p - q z 


‘s CA 


+ q z d sincj) - Zq^d cos tp 

+ 2q g d5> sincj) - i s x CA 

- q d coscj) 
n s r 


( 11 . 21 ) 


Finally, the quantity q^ can be related to q g by 


q = q +0 


and 


( 11 . 22 ) 


q = q + <P 


(11.23) 
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The procedure now will be to introduce the expressions derived for 


a^ , a^ , and q into equations (II. 7), (II. 8), and II. 9) in order to 
x y p 

develop expressions for A , A , and and a differential equation for <f>. 

x y 

Equation (II. 7) becomes 


- ^ = “ + % (Z CA • d slrll|,> 

P . . 

+ 9 S (w - 2d sintj) - 2d<J> cost})) 

+ (+ d cos4> - 2d<f) sln<\) -d<f> cosfy) 
+ q* ("X CA - d cost})) 

+ <}>(- d sint}>) 


(11.24) 


Similarly, equation (II. 8) becomes 

w + q g (-x " d cost})) 

• • 

+ q (-u - 2d cost}) + 2dtj> sine})) 
s 

+ (-d sine}) -2dt}) cost}) + d<}) sint})) 

+ q* (~z CA + d sint})) 

+ <K-d cost})) (11.25) 

Equation (II. 9) becomes 

- M. - A d cost}) - A d sint}) = I q + I t}> (11.26) 

A z T x yys s yyp 

Equations (II. 1) , (II. 2), (II. 3), and (11.26) are a set of four simultaneous 

first order differential equations in the four variables u, w, q , and 0. 

s 

However, they are not in a form amenable to numerical integration and equations 




(11.24) and (11.25) must be used along with considerable algebraic manipu- 
lation to put them in such a form. In order to simplify this process, the 
following definitions are made. 


b l = V Z CA ' d Slr " ( ’ ) 

b^ ~ -m (d sin$) 

b~ = (q (w - 2d sin <p - 2d(p costp ) 
j s 

+ (d cos (J) -2d(?Gin(J) - 2d(p 2 cos<p) 

+ q*(- x CA - d cos<J))} m p (11.27) 

With these definitions, equation (11.24) becomes 


• • 


-A = m u + b n q + b„d) 4- b~ 
x p l n s 2 3 


(11.28) 


Substitution of this into equation (II. 1) results in 


T -mu- b,q - b„d> -b_=mu + mqw 
x p l^s 2 r 3 s s^s 


(11.29) 


Equation (11.29) can be rewritten as 


(m + m ) u + 0 w + b 1 q + b_0 = T - m q w - b_ (11.30) 

s p X S Z X S S J 


In a similar manner, the definitions 
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C 1 = X CA ” d cos ^ m p 

^2 = m p ^ -d cos< £) 

C_ = {(-u -2d cosi}) + 2d_4> sintfOq (11.31) 

j S 

+ (-d sine}) — 2dcj> cost}) + d(J> 2 sin<f>) 

+ q 2 (~z CA + d sint}))} m p 


lead to rewriting equation (11.25) in the form 


-A — m w + C.q + C_(|) + C_ 
z p Is 2 T 3 


(11.32) 


Substitution of equation (11.32) into equation (II. 2) results in 


• • 


T - m w - CLq - C 0 d) -C 0 =mw-mqu 
z p Is 2 T 3 s s^s 


or 


0 u + (m + m )w + CLq + C 0 <j) = T - C, + m q u 

S p i. S M Z j s s 


(11.33) 


Next, it is necessary to substitute equation (11.28) and (11.32) into 
equation (II. 3). However, prior to this step, let us assume that is of 
the form 

M a = K(<f> - <J> o ) + C$ (11.34) 

This will make equation (II. 3) and (II. 9) independent so that neither has to 

be used to solve for M. to substitute into the other. With the indicated 

A 

substitutions, equation (II. 3) becomes 
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M A " 2 CaV - Z CA b l“s ' Z CA b 2 * 


' z CA b 3 + X CAV + X CA C l q S 


+ X CA c 2 * + X CA C 3 + M T " V= 


Regrouping the terms in equation (11.35) leads to 


<z CaV u + a yyx + z CA b l - x CA C l )<l s 


+ ( - x CA m p )w + (z CA b 2 - x CA C 2 )<t ’ 


- tt], + M A + x CA C 3 - z CA b 3 


Substitution of equation (11.28) and (11.32) into equation (11.26) yields 


M. + (m d cos<J>)w + (C.d cos4>)q 
A p 1 n s 


+ (C 2 d cosc{))<|) + C^d cos<{) 


+ (nipd sinc}))u + (b^d sin<{))q s 


+ (b 2 d sin<{>)<J> + b^d sin<f> 


= I q + I <f> 

yyp s yyp Y 


Collection of terms in equation (11.37) results in 
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(-in d sind>)u + (-m d cosrf>)w 
P P 


+ (I - b,d sind) - C,d cos<f>)q 
yyp 1 1 r s 


+ (I - C_d cosd) -b„d sin <b) <f> = 

yyp 2 T 2 r 


- + C^d cosd> + b^d sine}) 


Equations (11.30), (11.33), (11.36), and (11.38) are now of the form 


B ll u + % w + B 13 q s + B W* = D 1 


B U u + B 22 W + B 23 q s + B 24^ ' B 2 


B 31 u + B 32 w + B 33 q s + B 34^ " ®3 


B 41 U + ®42 w + B 43 q s + B 44^ * D 4 


where 


B 11 


B 


12 


m + tn 
P s 

0 


B 13 = b ] 
B i4= b : 


T - m q w - b„ 
x s s 3 


B 31 = Z CA m P 
B 32 = - X CA m p 

B 33 I yys + Z CA b l X CA^ 
B 34 = Z CA b 2 “ X CA C 2 
D 3 = M T + M A + X CA C 3 - z 


B„. = 0 


21 

22 

23 

24 


(m + m ) 

P s 

c i 

C 2 

T 4- m q u 
z s s 


- C. 
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42 


B 43 = 

B 44 " 
D, = 


-m d sind) 

P Y 
-m d cosd) 

P 

I -b-d sind) -C..C 
yys 1 T 1 

I -c„d cosd) -b„c 
yyp 2 2 

-M^ + Cgd cosd> + b. 


(11.38) 


(11.39) 


CA b 3 


. COS d> 

sind) 
d sine}) 
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The differential equations of motion for the system can be written as 
(let 0=0)) 





• \ 

u 


f \ 

D 1 




• 

w 


D 2 

B 


- 

• 

q s 

' = 

D 3 


- 


• 

0) 

W 4 


D 4 

V. / 


(11.40) 


plus the equation 


0 = 0 ) 


(11.41) 


Using a linear systems solver or matrix inverse at each derivative evaluation, 
the vector 


u 


w 




can be formed. 


Given initial conditions 


\ 

U 


f 

0 

W 


0 

• 




» = « 

Of 

• 


0 

03 





ie ; 


o. 
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an initial value for 0, i.e., 0 , values for d, d, and d, m , m , I , 

o ps yys 

I p, the spring constant K, the geometric quantities r^ (x^ ant * Z QA^ 
r CT ^ X CT an< * Z Cp’ anc * t ^ ie t h rust components T x and T , then it is possible 
to carry out the required numerical integration. 

As the integration is carried out, A , A , and m can be evaluated at 

X z 

each point in time from equations (11.28), (11.32), and (11.34), respectively. 



III. TWO JOINT MODEL 


In the two joint model, single axis pin joints will be assumed at the 
shoulder and the wrist of the shuttle/arm/payload combination. 

The forces and moments acting on the shuttle itself are shown below 
in Figure III.l. 



Figure III.l: 


Shuttle Free-Body Diagram 


where 


T = 

T i 

+ 

r1° + 

T k 

A = 

X 

A I 

+ 

'y 

z 

A k 

ii 

o 

i> 

X 

u i 

+ 

o 

+ 

z 

w k 

r CA ! 

= X CA 1 

+ 

0 

7^ + 

2 ca E 


_ /tO _ rO 




The quantities marked through are zero because the model is two-dimensional. 
The xyz coordinate system is centered at the shuttle center of mass, C, 
and rotates with the shuttle. 

The force and moment equations for the shuttle in body fixed axes are, 
as before [see equations (II. 1 thru 11,3) 3 . 


£F = T 

+ A 

= m (u + q w) 

(III.l) 

X X 

X 

s s 


ZY = T 

+ A 

• 

= in (w - q u) 

(III. 2) 

y y 

y 

s s 


£M_ = M. 

+ 

• 

+ z„ . A - . A = I q 

(III. 3) 

C A 

T 

CA x CA z yys s 



where all variables are defined as they were in the analysis of the single 
joint model. 

Equations (III.l) through (III. 3) can be written as differential equations 

• • • 

for u, w, and q with unknowns A , A , and M appearing on the right hand 

S X Z A 

sides. Thus, 

T A 


• 

u = -q w 


(III. 4) 

n s 

m m 

s s 

T A 

w = q u 

+ — + ~ 

(III. 5) 


m m 

s s 

• 1 

q s = I 

{K T + z CA A x - x CA A z + M A } 

(III. 6) 


yys 


As before, the torque exerted by the arm on the shuttle at point A will be 


modeled as 
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In order to find the forces and A , it is necessary to consider the 
arm and then the payload. 

The arm will be assumed to be massless. Thus, its effect is to transmit 
and modify torques. The free-body diagram of the arm is shown in Figure III. 2. 



Figure III. 2: Arm Free-Body Diagram 


The equations of motion for the arm are 


EF « -A + B = 0 

X XX 


EF = -A + B = 0 
y z z 


EM a = -M^ + Mg - B z d cos* -B^ d sin* = 0 


From equations (III. 8) through (III. 10), we get 



H = M. + A d cos* + A d sin* 

D A Z X 


(111. 8) 

(111. 9) 

(111. 10) 

(111. 11) 

(111. 12) 


(III. 13) 
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However, as was the case with the shoulder joint, a form will subsequently 
be specified for (in terms of an angle yet to be defined) . 

The next step is to examine the free-body diagram of the payload shown 

below 



Figure III. 3: Payload Free-Body Diagram 


The equations of motion for the payload are given by 


where 


£F = -B = m a^ = -A 

X X P D X 

* V 


EF = -B = m a_ = -A 
Z Z p D z 

* <7 


ES D - -V + ? DB * ( - 5) ■ Sjp V 


r DB r BD 

J db * ( - 5) 


+£sinijj i + £cosi|j k 


r BD * B 


r x A 
BD 


r DB X ( ' B) 


= (-£ sinipi - JicosiJj k) x (A^i + A_k) 
= (+ A Z simp - A Z cosiJj)j 

Z X 


(111. 14) 

(111. 15) 

(111. 16) 


x 
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Thus, equation (III. 16) can be written as 


% - ’“a + A z £sini|J - A x Icosip = ^yyp^p 


(III. 17) 


In equations (III. 14) and (III, 15), a^ and are expressed in a 

x z 

coordinate system which rotates with the shuttle (at an angular rate given 
by fi = w a = q g j). 

As before [see equation (11.13)] 


*D - 3 C + ®D , + 5 * (Si * ; CD> 
rel 


+ 2fi x v , + S2 x r__ 

rel CD 


(III. 18) 


In order to generate a complete expression for a^, it is necessary to first 
write an explicit form for r . Thus, 

vD 

r = r + r + r 
CD CA AB BD 

or 

r CD = (x CA 4- d cos<j) - £sini|j) i + (z^ -d sintj) - £cosi|Ok (III. 19) 

As in the single joint model, the length of the arm will be modeled 


N 


as 


d(t) 


d + d t + d t z /2 
o o o 


(III. 20) 


where d^, d Q , and d are input parameters. Now let us proceed to develop each 
of the terms of the right-hand side of equation (III. 18) individually. 
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From equations (III.l) and (III. 2), a £ is given by 

a =(u + q w)i + (w - q u)k (III. 21) 

C S S' 

Next, the quantities Vp and a^ will be developed from r^. Recall 

rel rel 

that r CD is given by 


r = (x + d cost}) “ Asin^)i + (z rA - d sint} - £cosl|))k (III. 22) 

CD CA CA 

where x , z , and Z are constants. To find V and a , i and k are 

vA wA " i i/ i 

rel rel 

treated as constant vectors also. Thus, 


V n = (d cost}) - d0sin<£ - Zipcosip)i + (-d sin<j> - d<}>cos(}) + Zip sin^)k (III. 23) 

rel 


and 

a_ = (d cost} - 2dt})sin<}> - dt}>sin<f) - d<}> 2 cos<p - Zipcosip + Zip sini|j)i 

rel (III .24) 

+ (-d sin({) - 2d0cos(j) - dc^coscf) + dt}> 2 sint}) + Zipsinip + Zip cosi|0k 


The quantity ft x (ft x r^) is given by 

q I x [<1 J x { (x r A + d cost} - £sini}))i + (z pA - d sint}> - £cost|))k}] 

S S CA CA 


Thus, 


S x (n x r CD ) =(-qgX CA - q 2 d cost}) + q 2 £sini}))i 
+ (-q 2 z + q Z d sine}) + q 2 £cos^)k 

S CA S S 


(III. 25) 


The term 2ft x is generated from 


2ft x V - = 2q j x {(d cost}) - d^sint}) - Zipcosip)i 
it ex s 

+(-d sint}) - dt}>cos<}> + £$sin^)}k 


and 



+ (-2q d cos<|> +2q d^sintj) + 2q JhjjeostjOk 
s s s 

Finally, the £2 x r^ term is generated by 

ft x r CD = q g J x {(x CA + d cos^) - £sintjj) i 

+ (z^ A “ d sin$ - £cosif>)k 


(III. 26) 


or 

Z. — • • • — 

x r CD a (q s Z CA “ <l s d sin<J) lcos ^ i 

+ (-q x p . -q d cos(j> 4- q £sini|j)k (III. 27) 

S v»A S S 

Combination of Equations (III. 21), (III. 24), (III. 25), (III. 26), and 
(III. 27) to form results in 



u + q g w + d cos<f> -2d<f)sin<f)~d(f)sin<f) 

- d4> 2 cos4> - £$cosljJ + £^ 2 sin ijj 

- q^x CA -q 2 d cos^ + q 2 S-sin^ 

• • • 

- 2q d sincj) -2q $cos$ + 2q £\jjsin^ 

s s s 

• • '• 

+ q z r . -q d sincf) -q £cosi|> 

S l»A S S 


• • • 

= w — q g U - d sin$ - 2d(pcos(p - dtficos^ 

+ d<{) 2 sin(J) + £ijjsin ip + £i|i 2 cosi[) 

- q 2 z p . + q 2 d sin(j) + q 2 £cost|J 

S LA S S 

- 2q o d cos 6 + 2q d(j) sin<f> + 2q £ij;cosij> 

s s s 

- q s x CA “ % d cos< P + q s ^sini{j 


(III. 28) 


(III. 29) 


Substituting equations (III. 28) and (III. 29) into equations (III. 14) and 


(III. 15) and collecting terms, we get 
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-A = m u + 0 w 
x p 


and 


-A = 


. • 

+ m p( Z Q A ” d sint}) - J£cosiJ;)q s 

+ m^(-d sintj))t}> + m^(~Zcos\p)\p 
+ m [q 2 (-x - d coscf) + £sim p) 

p S L*A 
«• 

+ (d cost}) -2d4>sin(j) -d(j) Z cos(p + £tj) 2 sim [>) 
+ q g (w -2d sin<|) -2dt}icos<}> + 2Z\psimp)] 


0 u + m w 
P 

+; m (-x^ -d cost}) + £sinljj)q s 

+ m^C-d cos<J>)4> + in (£sinljj)t}j 

+ m [q 2 (-z + d sint}) + £cosi|;) 

p S LA 

+ (-d sint}) -2dcf)cos(}> 4- d<}) 2 sin<J) + Zip z cosi(j) 
+ q g (-u -2d cost}) + 2dt}>sin({) + 2£tjJcost}i) ] 


With the definitions below, the expressions for A and A can be 

X z 


simplified. 


/ 

b^ = (z CA ~ d sint}) - £cost}))mp 

b 2 - nip(-d sint})) 

* * • ' 

bg = nip[q g (w -2d sint}) -2dc})cost{) + 2£i})sim})) 

+ (d cost}) — 2dcf)sin4> -dtj) 2 cost}) + £i}> 2 sintfj) 
+ ^g(“ x CA " d cost}) + £sini})) ] 
b^ - mp(-£cost})) 


(III. 30) 


(III. 31) 
greatly 


(III. 32) 
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= ( _X CA -d cos4> + X,sini|j)mp 

c„ = m (-d cos(j)) 

/ P 

• • • 

= m [q g (-u -2d coscp + 2d({>sin(j> +2 £^cosi|j) 

• • »o *5\ 

4- (-d sin<}> -2dcj)coscJ) 4- d(j) sine}) + to|j costy) 


+ qjj(-z^ + d sin<{) + Acosip) ] 


= m (&sini{j) 


(III. 33) 


The resulting expressions for A x and k^ are 


-A = m u 4- b,q 4- b„d) 4- b , xb 4- b_ 
x p Is 2 4 3 

-A z - m p w + + c 2 * + c 4 i + c 3 


(111. 34) 

(111. 35) 


Rewriting equations (III.l) and (III. 2), we have 

mu-A = T - m q w (III. 36) 

s xx s n s 

mw-A=T4-mqu (III. 37) 

s z z s n s 

Combination of equations (III. 34) and (III. 36), (III. 35) and (III. 37), 
results in 


^ (m s + V" + Vs + V + V ‘ T x * m s V ' b 3 

< m s + V" + c l% + c 2 + V ‘ T z m W ‘ c 3 
Equations (III. 3) can be rewritten as 


(111. 38) 

(111. 39) 


I q 4- x A - z_.A = M 4- M. 
yys^s CA z CA x T A 


(III. 40) 


Combination of equations (III. 34), (III. 35), and (III. 40) results in 
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I q - x_.m w - x-,.c.q - x„.c 0 <}> 

yys M s CA p CA l*s CA T 


Vi* - X CA C 3 + Z CA m p U + Z CA b l q s 


+z CA b 2 ‘ t ’ + z CA b 4 ^ + z 0 A b 3 ' “t + M A 


which can be rewritten as 


(z CA m p )u + ( ' x CaV " 


+ (I + Z A b- - x r ,.C 1 )q 
yys CA 1 CA 1 n s 


+(z CA b 2 ' X CA C 2 )( ^ 


+<Z CA b A ' X CA°A^ 


III. 41 


“t + M A “ Z CA b 3 + X CA C 3 

Combination oi equations (III. 34 , III. 35 ) and (III. 17 ) results in 


I q - m ftcosiiu - b^Jlcosi^q 

yyp p p 1 s 

- b 2 &cos# - b^Zcosipty - b 3 ilcosT|) 
+ mp&sinijw + c^£sinlpq s 
+ c^&sin# + c^£sin^ + c^^sin^ 


( 111 . 42 ) 


Finally, the quantity can be related to by 


q p = q s + <5 + (i> - }) 

q = q + jp 

q p = q s + $ (III. 43) 

Combination of the last of equations (III. 43) with (III. 42) and rearranging 
terms results in 

(- 111 ^ £cos^)u + (m £sin40w 

+ (lyyp “ &cos\jJ + c x £sini[i)q s 

+ (-^2 Icosifj + JlsiniJj)^ 

+ ( lyyp - b^ Zcosip + Isin^)^ 

= - + b^ &cost(j -c^ ilsiniji (III. 44) 

Equations (III. 38), (III. 39), (III. 41), and (III. 44) are four equations o'f 
the form 

B il“ + B i2^ + B i3 % + B i4^ + B i5^ = D i < i=1 » 4 ) (III. 45) 

This is a system of four equations with five unknowns. 

Equation (III. 13) can be used, along with (III. 34) and (III. 35) to 

obtain the needed additional relation. Equation (III. 13) can be rewritten as 

shown below when A and A are eliminated. 

z x 
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which can be rewritten as 


(tn^d sin<}Ou + (m^d cos$)w 

+ (b T d sintj) + c d costfOq 
X X s 

• • 

+ ^d sin<{> + c^d cos<t>)<J> 

• • 

+ (b^d sin^i + c^d cosc}))^ 

= M. - WL - b„d sincf) - c„d cos<j> 
A iS X J 


(III. 46) 


Equation (III. 45) forms the fifth equation needed. Now if we let 

* * “4, 

* - ; v 

and add the two differential equations 

5 = ^ 

* - 


we will have seven first-order ordinary differential equations in seven unknowns. 

The only step left is to specify the form of must depend on the 

relative angular displacement of the payload with respect to the arm (measured 
from its initial relative angular displacement, and on the angular rate of the 
payload with respect to the arm. The geometry for this is shown in Figure III. 4. 



Figure III. 4: Payload Angle Geometry 


39 • 


The required angle is a where 
a = 4> - iJj 


Also, 


and 


a = <J> - ib 

o o Y o 


* • • 

a = ij> - tJj 


The torque. Mg, on the arm, is negative when a > a Q and a > 0. Thus, 


M b = -K b (a - a 0 ) -C B a 


or 


Mg = -Kg (4, - ip - <j, Q + iP Q ) -Cg(<J> - i|>) 


(III. 47) 


The resultant elements of the b and D matrices arc: 


B- - = (m + m ) 

11 s p 

E 12-° 

B 13 * b l 

B 14 ’ b 2 

B 15 ' \ 

D 1 ' T x - m s q s K ' b 3 


B A1 = -m^Jicosijj 

a ! 

P 


41 p 
B ^2 ~ m &sin\}i 


B /0 = I -b,£cosi}j + c..£sind/ 
43 yyp 1 1 

B^ = -b2&cosi|j + C2-?-sin^ 

B._ = I -b.£cosdj + c.AsiniJ; 
45 yyp 4 4 

= -Mg + b^Jlcos^ -CgJlsin\[) 


E 2l'° 

B 22 = (n s + V 

B 23 " C 1 
B 24 ' C 2 
B 25 ” C 4 

D 2 =T z + m S V C 3 


B_, = m d sind) 

51 p T 

B__ = m d cosd) 

52 p 

= b^d sin<J) + c^d cost}) 

B^ = b 2 d sintj) + C 2 d cost}) 

B__ = b.d sind) + c.d cosd) 

55 4 T 4 

Dj. = M a - Mg - b^d sin<{) - C^d cost}) 
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IV. NUMERICAL RESULTS 


The one joint model developed in Section II and the two joint model 
developed in Section III have been compared for a series of thruster firing 
inodes. The goal of this comparison was to determine the similarities and 
differences in the dynamics as predicted by the two models. An analysis made 
with the programs written to check the algebraric signs and the initial 
magnitudes (see Appendix A) indicated that for the loaded RMS extended vertically 
above the shuttle, Thruster 9 (the center RCS thruster firing upward on the 
shuttle nose) produced the largest forces and bending moments at the shoulder 
joint. All thrusters were checked in this determination. 

Three comparison cases are presented here as an example of the data 
obtainable from the programs. The three cases assume that Thruster 9 is fired 
for 0.1 seconds and investigates the system oscillation during the subsequent 
100 seconds. The parameters for each of the three cases are listed in Table 1. 

The numerical results were obtained by numerically integrating equations 
(11.40) and (11.41) for the single joint model and equations (III. 45) and (III. 46) 
for the two joint model using the parameters listed in Table 1. The equations 
were integrated with a 4th-order Runge-Kutta integrator. The simulation assumes 
that the RMS is manipulating a cylindrical payload with a mass of slugs 

The numerical results obtained in this study are shown in Figures IV. 1 through 
IV. 6. Two separate cases were considered with the two-joint model. The cases 
differed in the length of the wrist joint. In Case II, the wrist joint was 
zero while in Case III the wrist joint was 7.5 ft. 



Table 1. PARAMETERS FOR NUMERICAL SIMULATION 


Program 

m 

P 

m 

s 

d 

• 

d 

d 

X CA 

Z CA 

C A 

k a 

I YYS 

I YYP 


Case 1 
Single Joint 


Case 2 
Two Joint 


w 

o 

^so 

00 =00 


(j)0 


0) 


1^0 

C B 


2018.5 slugs 
4689.4 slugs 
50 ft 
0 
0 

36.33 ft 
-6.083 ft 

1.22 x 10 5 lb. ft . sec/rad 

1.0 x 10 4 lb. ft/radian 
5699753. slug. ft 2 
625805. slug. ft 2 

0.0 lbs 

870.0 lbs 
0 

0 

0 

0 

90° 

0 

0 

0 


Same as Case 1 


1.22 x 10 5 lb. ft. sec/rad 
1.0 x io 4 lb. ft/rad 


PLOT SYMBOL 


0 


Case 3 
Two Joint 

Same as Case 1 


7.5 ft 
0 
0 

Same as Case 2 
Same as Case 2 

<$> 
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The data for Case 1 do not contain the wrist joint rate, joint angle, 
or moment at the wrist because these variables did not appear in the single 
joint model. These variables do, however, appear in Cases 2 and 3 which 
were generated using the two joint model. The output variables (all plotted 
versus time) for the three cases are as listed below. 


SH VX 


SH VZ 


QS 


PHIDOT 


PHI 

ALPHA DOT 


ALPHA 


ARMFX (=A X ) 


This is the x component (in body axes) of the 
velocity of the center of mass of the shuttle in 
This is the z component of the velocity of the shuttle 
center of mass in ft/s. 

This is the pitch rate of the shuttle in degrees 
per second. 

This is the angular rate of the shoulder joint 
in degrees per second. 

This is the shoulder angle in degrees 

This is the wrist joint rate (rate of rotation of the 

payload relative to the arm) in degrees per second. 

This variable is not integrated but is created after 
integration by differencing angular rates 
This is the wrist joint angle in degrees. This 
variable is the angle between the arm and the payload 
x axis. 

This is the x component of the force in pounds at 
the shoulder exerted on the shuttle by the RMS (The 
x component of force exerted on the payload by the wrist 
is -AX for an RMS which is modeled as massless). 
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ARM FZ(=A^) This is the Z component of the force in pounds at 
the shoulder exerted on the shuttle by the RMS. 

AM(=M.) This is the shoulder moment exerted by the RMS on the 

A 

shuttle in ft. pounds. 

BM(=Mg) This is the wrist moment exerted on the RMS by the 

payload in ft. pounds. 

Figures IV. 1 through IV. 6 give 100-second time histories of all output 
variables for the three cases described previously. Recall that case 1 does 
not output the variables a (ALPHADOT), a (ALPHA), and BM. 











PHI (degrees /sec) PHIDOT (degrees /sec ) 

89.79 89.89 89.99 90.09 90.19 <-“ 0.04 - 0.02 0.00 0.02 0.04 





ALPHA (degrees) ALPHADOT (degrees /sec ) 
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CVI 


FIGURE IV. 4 



TIME (SEC) 


i 

co 



TIME (SEC) 


FIGURE IV. 5 

O 



TIME (SEC) 







From the figures, it is apparent that the dynamic behavior of Case 1 
(symbol A') matches that of Case 2 (symbol 0 ) quite well. Both amplitudes 
and periods of oscillations match very well. This match is significant 
because the case 1 curves (A) were produced with the single joint model 
while the Case 2 curves (□) were produced with the two joint model. 

The similarities in dynamic behavior between Cases 1 and 2 can be at- 
tributed to the following factors: 

1. The payload center of mass is the same distance from 
the shoulder (50 ft) in both cases. 

2. The shoulder spring constant and damping constant are 
the same in both cases. 

3. The excitation (Thruster 9) is the same in both cases. 

4. The initial conditions are the same. 

The differences between Cases 1 and 2 are as follows. 

1. The wrist in Case 1 is rigid while there is a joint with 
a spring constant and dampoing at the wrist in Case 2. 

2. No wrist joint torque is obtained in Case 1 (it could be 

obtained by a rigid-body analysis but would be far too 
large), while the torque is obtained in Case 2. 

The results for Case 3 (symbol <£> ) do not match the results of Cases. 1 

and 2 mainly because the payload, c.g., in this case is 7.5 feet away from 

the wrist joint (the more realistic case). In Cases 1 and 2 the payload is 

assumed to have its c.g. at the wrist joint (in all cases the payload has a 
non-zero moment of inertia about its center of mass). 
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The period of oscillation in Case 3 is longer than in Cases 1 and 2 
(31 seconds as opposed to 27 seconds), largely because of the additional 
7.5 feet between the payload center of mass and the shoulder joint (the 
pendulum is longer). 

The magnitude of the oscillations in the x component of velocity of 
the shuttle’s center of mass is about the same in all three cases whereas 
the z velocity component oscillation is slightly larger for Case 3. The 
amplitude of oscillation in shuttle pitch rate is also larger for Case 3 
than for Cases 1 and 2. 

The angle <J> between the arm and the orbiter oscillates with about 
the same amplitude and rate variations in all three cases. However, integration 
over a longer time might reveal an energy exchange between the wrist and 
shoulder joints for the two joint cases (Cases 2 and 3). This should be 
looked into when the models are compared fully. 

The angle a and its rate a are defined only for Cases 2 and 3. 

After an initial jump of 0.05 radians during the thruster firing, the angle 

a oscillates over a range of about 0.1 degree in both cases. The differences 
in geometry between the two cases causes the curves for a and ct to differ 
between the cases. 

The crucial values in the analysis are the magnitudes of forces 
s 'and moments at the shouler (and wrist) joint(s). In all three cases, the 
force component in the z direction is the important force variable. Although 

the z component of force oscillates after the thruster is cut off in the same 

way that the x component oscillates, the oscillation cannot be seen on the 
plot because of the very high peak values of the z force which occur for all 
cases during the thruster firing (during the initial 0.1 seconds). The force 
peaks were in the range of 550 lbs to 600 lbs in all three cases. Due to the 
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geometry, this force is an axial force on the arm. This force acts at both 
joints (shoulder and wrist). For the values of system parameters chosen, 
the peak values of force at the joints occur during or at the end of the 

thruster firing after which the joint forces become very small. 

The moments experienced at the shoulder and wrist joints peak 90° out 
of phase with the forces at the joints. The plots of AM (all three cases) 

and BM (two cases) show how the moments vary as a function of time. The 

maximum torques exerted at the shoulder are larger than those at the wrist 
for the values of system parameters chosen for the examples. This is because 
the mass (and moment of inertia) of the shuttle are larger than those of the 
payload. 

Examination of the curve for the wrist moment (BM) indicates clearly 
the effect of payload rotation about the wrist joint. Case 3 with the payload 
c.g. offset from the joint gives higher values of wrist moment as would be 
expected. Also, the predominant period of oscillation is again longer for 
Case 3 than for Case 2. 
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V. SUMMARY , CONCLUSIONS AND RECOMMENDATIONS 

Both single joint and two joint models have been developed for 
two-dimensional analysis of shuttle RMS joint* loads. A limited comparison 
of the two models has been made which indicates the following. 

1. In many situations, either the single joint model or the two 
joint model can be used to predict the approximate maximum 
magnitudes of shoulder joint loads (forces and moments). 

2. The single joint model can be used to predict maximum wrist 
joint forces, but the two joint model must be used if 
information concerning the wrist joint moment is desired. 

3. Both models give information on shuttle oscillations induced by 
the RMS /payload combination. The time histories of the oscil- 
lations at each joint are also available. 

Both the one joint and two joint models are limited. The major 
limitations are as follows. 

1. Both models are two-dimensional while the RMS/shuttle/payload 
system is a three-dimensional dynamical system. 

2. The number of joints in the models is one or two while the 
number of joints on the actual RMS is actually six. 

Although the two limitations listed above are major, the models developed 

should be quite useful because of the joint sequence on the RMS (three of 

the joints are collinear and occur in sequence - one at the shoulder, one 

at the elbow, and one at the wrist) . 

\ 

The procedure by which the models were developed is readily extendable 
to a three joint (elbow joint added) planar model. Such a model is the next 
logical step in building up to a simplified dynamic model of the three- 
dimensional shuttle/RMS/payload system which might be usable in a simulator 
(a full-blown dynamic model will have difficulties operating in real time, 
and real time is a must for training simulators). 



It is recommended that the following items be studied further. 

1. The three joint planar model should be developed. 

2. A two joint, three-dimensional model should be developed 
(shoulder yaw and shoulder pitch) . 

3. A four joint, three-dimensional model should be developed as 
an extension of the model of (2) above, and then a three- 
dimensional six joint model should be developed. The philosophy 
of each model should be to keep it as simple as possible while 
performing the required tasks. 

4. The models completed should be compared, and versions selected 
for use as needed in: 

(a) training simulators, 

(b) mission planning and sequencing programs, and 

(c) engineering simulators. 
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APPENDIX A 



ANALYSIS TO CHECK SIGNS OF A , A , etc. at t = 0 

x’ z 


SINGLE JOINT MODEL 

All rates will initially be zero, and quantities will be determined only at 
the instant when the thruster is turned on. Equations (II. 1) through (II. 3) 
are 

(11. 1) 

(11. 2) 


T +A = m (u 4- q w) 
x x s s 


T + A = m(w-qu) 
z z s n s 

m a + z ca a x - X CA A 2 + < XI - 3 > 


At t = 0, q g = W = u = 0. Also, (j) = j ,• and = 0. Thus, from equation (11.34), 

• • 

i.e., assume also that d = d = 0. 


M a = K(<p - <P Q ) 4- C(p 


(11.34) 


it can be seen that = 0 (at t = 0) . Thus, equations (II. 1) through (II. 3) 
become 


T 4- A = m u 
X X s 

T 4- A = m w 
z z s 

z . A - x A + M - I q 
CA x CA z T yys s 


f (A.l) 
(A. 2) 
(A. 3) 


From equations (11.28) and (11.32), we have 


-A = m u 4- b.q + b„d> 4- b„ 
x p l^s 2 T 3 

• • 

• * 

-A = m w 4- c,q 4- c„d> 4- c„ 
z p l n s 2 T 3 


(11.28) 

(11.32) 


where, from equations (11.27) and (11.31) give b^, b 2 > b^ and c^, c^, c 
respectively. At t = 0, these values are 


A-l 



A-2 


b 

b 

b 

c 

c 

c 


1 

2 

3 

1 

2 

3 


= m p (z CA " d Sin(i>) 

= -m d sind) 

P 

= 0 

- ( — x . ~ d cosd))m 
CA p 

= -m (d cosd)) 

P 

= 0 


Thus, (11.28) and (11.32) can be rewritten as 

-A x = m p u + b 1 q s + b 2 <f> (A. 4) 

-A z = m p w + + e 2 4> (A.5) 

Equation (11.38) can be written at t = 0 as 

• * 

(-m d sind>)u + (-m d cosd>)w 
P P 

+ (I -b,d sind) - e,d cosd))q 
yyp 1 1 T 

+ (I -c„d cost}) - b_d sind))d) = 0 (A. 6) 

N yyp 2 2 r T -i. 

Equations (A.l) through (A. 6) form a system of six linear equations in six 
unknowns . 

The equations, in matrix form, can be written as shown on the following 


page. 






TOO JOINT MODEL 


Equations (III.l) through (III. 3), with the assumptions q = u = w = 0, M. = 0 

S A 

gives (also d = d = 0) (tj) = ij; = 0, <}> = <j> ,t|> = 


m u - A * T 

s x x 

m w - A = T 

s z z 


I q + x_.A - z_.A = M_ 
yys s CA z CA x T 


Equation (III. 13) given = 0 at t = 0, gives 
A d sintj) + A d cost}) = 0 

X Z 

Equation (III. 17) can be written as 

I q + I i) + A £cosii> - A £sinU> = 0 
yyp s yyp x x 


(A. 7) 
(A. 8) 
(A. 9) 


(A. 10) 


(A. 11) 


Equations (III. 34) and (III. 35) can be written as 


m u + b.q + b_4> + b.ljJ + A 

p l^s 2 4 > 


-A 


9 •• •* 
m w + c n q + c„cb + c.\b + A = - s?l 
p l^s 2 T 4 Y z / 3 


(A. 12) 
(A. 13) 


Where, from equations (III. 32), (III. 33), and the assumptions, we have 



c 

c 

c 

c 


1 

2 

3 

4 


= m (z^ -d sine}) - £costy) 


= m (-d sind)) 
P 

= 0 


= m (-icosiJj) 

P 

= m (-x -d cose}) + £sin^) 

p wA 

= m (-d cosil)) 

P y 

= 0 

= m (Jlsind;) 

P 


NOTE: bp b 2> b^ and Cp c 2 , 
are formed as in TW0JNT. 




